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Introduction 


The goal is to softly land high mass vehicles (10s 
of metric tons) on Mars 

• Supersonic Retropropulsion (SRP) is a 
potential method of deceleration 

• Current method of supersonic parachutes 
does not scale past ~1 metric ton 

• CFD is of increasing importance since flight 
and experimental data at these conditions is 
difficult to obtain 

• CFD must first be validated at these 
conditions 
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Cianciolo, A. D., et al, 
"Overview Of The NASA 
Entry, Descent and 
Landing Systems 
Analysis Exploration 
Feed-Forward Study," 
International Planetary 
Probe Workshop, 2011. 



Fig. 3. Exploration Feed Forward Concepts 


• The EDL SA Team identified SRP as the only credible method of supersonic 
deceleration for Exploration Class (40 metric tons) vehicles entering Mars 
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SRP Flow Structure 


Supersonic Flow 


CFD of Daso et al (AIAA 2007-1423) 
Sonic nozzle 
M s = 3.48, C T = 0.4 


Bow Shock 


Sonic Nozzle 


Model and Sting 


Time = 0.000000 seconds 


SRP Flow Structure 


Supersonic Flow 


CFD of Daso et al (AIAA 2007-1423) 
Sonic nozzle 
M , = 3.48, C T = 0.4 



Opposing Jet 


Time = 0.000000 seconds 


SRP Flow Structure 


CFD of Daso et al (AIAA 2007-1423) 
Sonic nozzle 
M , = 3.48, C T = 0.4 



Time = 0.001035 seconds 


SRP Flow Structure 


Supersonic Flow 



CFD of Daso et al (AIAA 2007-1423) 
Sonic nozzle 
= 3.48, C T = 0.4 


Triple Point (2D) o 

Annular Rin? 

m m I I I I I V* " 


Recirculation 


Stagnation Region 


Time = 0.0000000 seconds 


CFD Validation Approach 


Employ multiple solvers to the same SRP problems 

• DPLR (Kerry Zarchi, ARC) 

• FUN3D (Bil Kleb, LRC) 

• OVERFLOW (Guy Schauerhamer, Jacobs/JSC) 

Compare results between codes and with historical tunnel data 

• Qualitative: Shock structure and standoff distance, unsteady behavior 

• Quantitative: Surface pressure, forces and moments 
Perform CFD-validation wind tunnel tests of SRP 

• Complete run conditions, quantified tunnel uncertainties 

• Higher thrust coefficients to better match flight requirements 


DPLR 

FUN3D 

OVERFLOW 



Verification 


Code-to-Code 

Comparisons 


Historical Tunnel Tests 


CFD Validation 
Tunnel Tests 
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Validation of 


Code-to-Test 

Comparisons 



CFD for SRP 


9 




Wind Tunnel Model 


• Air as freestream and jet gas. 

• 4 removable nozzle plugs. 

• 167 pressure taps including 7 
high frequency pressure 
transducers. 

• High speed Schlieren video 
(10 kfps). 

• Same model used in both tests. 



Forebody Instrumentation Aftbodv Instrumentation 

o ESI 1 Ports 



180 


• Kuliles 


x/L 



9 

(deg) 
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Wind Tunnel Tests 



Langley 4'x4' Test 

• Mach 2.4, 3.5, 4.6 

• Cj's up to 3, a couple at 6 

Ames 9’x7’ Test 

• Mach 1.8, 2.4 

• Cj's up to 10 

• Liquefaction in plumes 


LRC 


ARC 





Langley Test 






OVERFLOW Best Practices for SRP 


• Best practices based from LRC UPWT Run 165: 1-nozzle, Mach 4.6, C T 2 

• Grid refinement and time step sensitivity studies 

• Grids between 80 and 90 million points 

• Time steps <= 1.71e-06 seconds 

• 5 Newton subiterations 

• HLLE++ numerical flux function with Van Albada limiter for spatial terms 

• Symmetric Successive Over Relaxation (SSOR) algorithm with Newton dual time stepping for temporal terms. 

• Used Direct Eddy Simulation (DES) turbulence modeling with Menter's Shear-Stress Transport (SST) as the Reynolds 

Averaged Navier-Stokes (RANS) submodel. 


• For SST, used strain-based production term employing Wilcox's realizability constraint 
• All jet-on cases were solved time-accurately. 




Structured Overset Grid System 





US3D FUN3D TEST 


LRC Run 165 Qualitative Comparison 
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1-nozzle, Mach 4.6, C T 2. 
Periodic unsteadiness. 

Oscillating pressure 
wave increases average 
pressure on model face. 

Movies not synced in 
time. 


US3D solution by Emre 
Sozer (ARC) 

Cart3D solution by Noel 
Bahktian (Stanford) and 
Michael Aftosmis (ARC) 


cc 

< 

o 


15 





LRC Run 165 Comparisons 


0 . 06 - 


DPLR SST 
FUN3D SA-DES 
OVERFLOW SST-DES 
US3D SA-DES 
Test 


0 . 04 - 


0.02 - 



- 0.02 - 



Capturing periodic 
oscillation in the triple 
point increased average 
pressure on the face. 

All codes fall within 
tunnel uncertainty. 



Turbulence Modeling 


Unsteady behavior was influenced by 
turbulence modeling. 

• For SRP, limiting eddy viscosity 
produced more realistic behavior 

• Cart3D is inviscid 

• Figure is ratio of eddy to laminar 
viscosity for different turbulence 
models, simulated with 
OVERFLOW and FUN3D. 

• For SRP simulations: 

• DPLR used SST-V. 

• FUN3D used SA-DES 

• OVERFLOW used SST-DES 
and SST-RC 



(a) SST: Me nter r s baseline it-uj SST model. 4 ®’ 


(b) SST-CC: SST with Suzen and HofFmanrv& compressibility 
correction. 47 




(c) SST-V: SST with vocticity-based production . 30 


(tl) SST-RC: SST with WMcck's realizability constraint . 41 



(e) SST-DES: SST-based DES . 42 


(f) SA-DES: SA-based DES . 34 


Figure- 20: The ratio of turbulent eddy viscosity to laminar {bulk} viscosity for various turbulence models. 

Note: No attempt was made to capture these instantaneous snapshots near the same point of the quasi-period ic cycle. 
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LRC Sting Sensitivity Study, 3-nozzle 



0 0.5 1 1.5 2 

x/D 


OVERFLOW FUN3D DPLR TEST 


LRC Run 165 Comparisons 








LRC Run 165 Comparisons 



Pa, total ~ ^A,aero + ^l 







20 


OVERFLOW FUN3D DPLR TEST 


LRC Run 262: 3-nozzle, Mach 4.6, C T =3 
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LRC Run 262: 3-Nozzle, C T =3, a=12° 



- 0.05 

- 0.1 



Bow shock 
shedding impacts 
model face and 
side shell 

Constructed 
shadowgraph of 
CFD solution 
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Run 262: 3-Nozzle, 


aO° 





C t = 3, <t>=0° 





DPLR steadier than 
FUN3D and 
OVERFLOW (SST- 
V vs DES) 


Large scatter in 
neighboring 
pressure ports on 
the mode 

windward side shell 


FUN3D 

overpredicts C P on 
model face for a=0° 
and a=12° 


Deviation at nose 
implies jet-to-jet 
interactions 
predicted differently 
between codes 


Deviation at 
shoulder implies 
differences in shock 
shedding impacting 
the model face 
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OVERFLOW FUN3D DPLR TEST 


Run 307 and 311: 4-Nozzle, C T =2, <j>=0°, 180° 


• Runs only differ in roll angle. 

• Short blunt behavior vs. large shock 
standoff. 
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Viscous Full Tunnel 


Triple Nozzle 
Run 267 

LaUPWT Section 2 
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Ames Test 




Z (diameters) 


Ames Sting Sensitivity Study, 1-nozzle 



Ames UPWT 
Single Nozzle 
Mach 2.4, C T 6, Re/ft 1M 


- 2-1012 

X (diameters) 


Time = 0.000044 seconds 




Run 223, (3 = 0°, 4°, 8°, and 12° 


Run 223 
1 -nozzle 
Mach 2.4, C T 4 


Run 223 
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Mach 2.4, C T 4 

3 4° 
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Run 141, (B = 0°, 4°, 8°, and 12° 


Run 141 
3-nozzle 
Mach 2.4, C T 6 


Run 141 
3-nozzle 
Mach 2.4, C T 6 

(3 4° 


Run 141 
3-nozzle 
Mach 2.4, C T 6 

P 8° 


Run 141 
3-nozzle 
Mach 2.4, C T 6 
P 12° 
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3-Nozzle, Mach 2.4, C T 10, Run 145 



KM 



Run 130, C T 6, (B = 0° 


• Comparisons with other codes 

• DPLR reached steady state, captures side shell C p well. 

• FUN3D and OVERFLOW simulations very similar 

• OVERFLOW captures C n and nose and near nozzle 

r 




Time = 0.000044 seconds 


c p 

■ 0.3 
0.15 
0 

- 0.15 
™ - 0.3 


NASA Ames 9x7 Supersonic Wind Tunnel 
Entry, Descent, and Landing Project 
Supersonic Retropropulsion Test 97-0234 
Scott.A.Berry@nasa.gov, Matthew.N.Rhode@nasa.gov 
August 2011 


Time 0.00 ms Recorded at 20,000 fps 







4-Nozzle, Mach 1.8 


Run 166 
4-nozzle 
Mach 1.8, C T 2 


Run 170 
4-nozzle 
Mach 1 .8, C T 6 


Run 172 
4-nozzle 
Mach 1 .8, C T 8 



■ 

■ 
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Thrust Dominance 


Aerodynamic forces are small when 
compared to thrust. 


Are aerodynamic forces negligible? 

• Need to know entry angles 
and vehicle design 

• Scalability 
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Computational Costs 

• From the Ames post-test runs 

• All cases run on Pleiades, either Nehalem or Westmere nodes 

• This is not a perfect comparison 

• DPLR is an average of three estimates 

• FUN3D is from a single case (not an average) 

• OVERFLOW is an average of six runs 

• Need time accurate runs for predictions 

• Would need to cut down computational costs for parametric studies or database 


generation. 

Solver 

CPU Hours 
per Case 

Iterations 
per Case 

Grid Points 

CPU seconds/ iteration/ 
grid point 

FUN3D 

28000 

39500 

42M 

6.1e-05 

DPLR 

44500 

106000 

53M 

2.9e-05 

OVERFLOW 

35039 

73500 

85M 

2.0e-05 
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Conclusions 

• Best agreement between WTT and CFD is for 1-nozzle cases and high thrust 3-nozzle 
cases. 

• Worst agreement was for low thrust 3-nozzle cases and high thrust 4-nozzle cases. 

• 3-nozzle cases more steady at C T s > 4. 

• 4-nozzle cases highly unsteady and chaotic at C T s >~ 4. 

• WTT averages spanning 2.5 seconds are probably not converged. 

• CFD averages spanning < 0.01 seconds are not too comparable to WTT data. 

• Large difference in data acquisition rates exists between codes and test. 

• Test rate was 10 or 30 Hz (0.1 or 0.033 seconds per reading) for 2.5 seconds for 
25 or 75 points per average. 

• CFD rate was between 190 and 400 points over a maximum of 0.017 seconds. 

• Frequencies captured by CFD not captured by test, and vice versa. 

• A "converged" average for CFD may be a completely different than the average 
that was obtained by the test data acquisition system. 

• Aero effects are small when compared to thrust. 

• Computational costs high for validation, could be much less for production. 
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Future Work 


• More steps need to be taken to better simulate 
Mars EDLSRP 

• NASA funding for SRP was discontinued in 
Fiscal Year 2012 

• Live rocket engine test including startup in 
SRP conditions 

• Sounding rocket test 

• CFD of flight conditions 

• Atmospheric and rocket 

• SRP is an enabling technology which still needs , 
development 

• Large-scale propulsion 

• Aero/aerothermal analysis 

• Vehicle design 

• GN&C 

• Additional funding avenues are being pursued 

• SpaceX and the USAF are researching SRP 
for Return To Launch Site capabilities 

• Masten Space is interested in returning 
rockets using SRP 
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